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Abstract 

Making use of the NJL model and the multiple reflection expansion approxima- 
tion, we study the phase transition of the finite size droplet with u and d quarks. 
We find that the dynamical masses of u, d quarks are different, and the chiral sym- 
metry can be restored at different critical radii for u, d quark. It provides a clue 
to understand the effective nucleon mass splitting in nuclear matter. Meanwhile, 
it shows that the maximal isospin chemical potential at zero temperature is much 
smaller than the mass of pion in free space. 



PACS numbers: 25.75.Nq, ll.10.Wx, 12.38.-t 
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The stability of quark matter has been investigated recently with various approaches, 
such as lattice QCD [H 0, MIT bag model 11 13 , NJL model EJ EDI [HI C2 • 
The strange quark matter containing u, d and s quarks are considered more stable than 
those only containing u, d quarks jU IS]- In the MIT bag model, hadrons consist of free 
(or only weakly interacting) quarks which are confined to finite region of space (in the 
"bag"). The confinement is not a dynamical result of the theory, but an input parameter 
with appropriate boundary conditions in the model. However, an important feature of 
MIT bag model is the absence of the dynamical chiral symmetry breaking. When the 
effect of the dynamical chiral symmetry breaking is taken into account, the strange quark 
matter is not absolutely stable Meanwhile, the transition from u, d quark to 

u, d, s quark by weak process is not favored, because when the chiral symmetry of u, d 
quark is restored, the symmetry of s quark is still broken ^Tj. It is then necessary to 
study the property of the system consisting of u and d quarks. 

So far all the foregoing analysis are either restricted to the quark matter that u, d 
quarks have the same chemical potential and the isospin is symmetric in the Lagrangian, 
or that of infinite volume. This means that, for the infinite quark matter, all the quantities 
related to u, d quarks in the matter are the same, in particular the dynamical quark mass 
and the quark condensates are equal for both flavors. However, there are many situations 
in the real world that the u, d quark numbers are not equal. According to the principle of 
statistical physics, the chemical potential should be different for u, d quark, respectively. 
For example, the neutron stars must be electrically neutral to a very high degree (even 
though there are electrons in neutron stars, the ratio is very small) J2]- Therefore if the 
neutron star consists of deconfined u, d quarks, the number of the d quark must be about 
two times of that of u quarks to ensure the electrical neutrality. On the other hand, the 
finite size effect including the contribution of the surface tension and the curvature, which 
has been shown to be very important to the properties of strangelets [UH], and of great 
interest in the study of nucleus, neutron stars, and heavy ion collisions JS], has not yet 
been taken into account for the isospin asymmetric quark matter [2^ ll6[lTT|IT%lll9|l2*Ul l21|. 
It is thus imperative to investigate the variation behavior of the isospin asymmetry and 
the chiral symmetry breaking in finite size quark matter (or quark droplet). Since the 
Nanbu-Jona-Lasinio (NJL) model [S] at quark level involves the information of 
chiral symmetry breaking and restoration, and of which the interaction among quarks 
is easy to deal with, in this paper we take the NJL model and the multiple reflection 
expansion (MRE) approximation [TUJ HU 122] to study the finite size effect on the 
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isospin asymmetric quark matter. For simplicity, we take only the u, d quark matter into 
account. 

In the two-flavor NJL model, the Lagrangian is written as 

£ = q{iYd^ - m + fir/ Q )q + L int , (1) 

where mo is the current quark mass. For the flavors u and d, one usually takes approxi- 
mation m u = m d = m . The interaction part is 

Hint = G[(qq) 2 + (#75 t*?) 2 ] . (2) 

In order to describe the isospin symmetry breaking, we introduce different chemical 
potential for u, d quark, respectively, 

V o n d J 

It is obvious that this breaks the isospin symmetry, and the breaking can be manifested 
by the isospin chemical potential /ij — /i u — fid, so different quark condensates (f> u = {uu) 
and (f) d = (dd) should be considered. The quark condensate is generally given by 

m = -*J (04^) , (4) 

where S(p) is the dressed quark propagator and the trace is on the Dirac and color spaces. 
After some derivations, the thermodynamic potential at finite temperature T and chemical 
potential /if (f — u,d) can be written as 

n(T,fi u ,fi d ,(fiu,(fid) = 2G((fi 2 u + (fi 2 d ) + ^M f (T,fif), (5) 

f=u,d 

with 

d 3 k 



/d k ( 1 
[E kJ + Tln[l + exp(--(E kJ - 

+Tm[l + exp(-i(£ fe , / + //,))]} , 



where E k j = \J M*j + k 2 , Mf is the constituent quark mass and can be given at the mean 
field level as 

M f = m - 2G<p f . (6) 
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For the u, d quark system in a spherical bubble with radius R (or referred as a quark 
droplet), the chemical potential is dependent on R, so we can write the Lagrangian as 



£ = g(i 7 ^9 - m + n(R)~f )q + £ 



inti 



(7) 



where p(R) = pp, the concrete form of p will be given later. Because the droplet (or 
bubble) has a finite surface and the curvature provides a pressure difference between the 
inner and the outer parts of the droplet, not only the quark matter in the droplet, but also 
the surface and the curvature contribute to the thermodynamical potential of the system 
[23] . To incorporate all these effects, an approach denoted as multiple reflection expansion 
approximation has been developed |BJ [TUl HH 122] • In the multiple reflection expansion 
approximation, for a droplet composing of quarks with flavor i, the thermodynamical 
quantities can be derived from a density of states in the form 



<m 
dk 
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k 2 V 



2tt 2 



+ fs(-)kS + f c (-)c + 



where V is the volume of the droplet, S = AttR 2 and C = 8irR are the area and the 
extrinsic curvature of the surface of the droplet, fs (^)> fc (^) are the contributions to 
the density of the states from the surface, the curvature of the droplet, respectively, and 
can be given explicitly as [HI El 



/ s (- 



arctan 



•57T 



7T 



(9) 



1 



12vr 2 



— — arctan ( — ) 



(10) 



3A; 

Then the density of states of the i flavor quarks can be given in the multiple reflection 
expansion approximation as 

12tt 2 

) - 



p = PMRE(k,m,R) 



Vffl./ 



kR 



(kR) 2 ^%n 



k 



(11) 



From this expression we know that, if we do not consider the contributions from the surface 
and the curvature, the Lagrangian in Eq. (J7J) has the general form used in studying infinite 
quark matter in the NJL model. 

By using the multiple reflection expansion approximation, we obtain then the density 
of states as [k 2 Pmre)I^ 2 and the thermodynamic potential can be given as 
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0(T , , rh rh \ ( M / ~ ^ 0M V- f PMRE^dk f 

, Pu, Pd, <Pu, <Pd) — }^ 2G ^ J 2tt2 1 hJ 

f=u,d - f=u,d 

+Tln[l+exp(-i(S Jfc>/ -/i / ))]+rin[l+exp(-i(S Jfc>/ +|i / ))]}.(12) 
In the zero-temperature limit T — > 

ftfT i i \ ( M /~ m o) 2 O/U f A E k,npMR E k 2 dk /*£ p MRE k?dk 

^ ^ = 2G " 2iVc 4 2~7 2 2iVc/iM 7o 

r A E k4 p M REk 2 dk f kd F PMREk 2 dk 



r &k,dPMREk dk [ k f pMREk dk 

L F — 2^ 2A Wo ' (13) 



-2iV c 



where p u = y(M u ) 2 + (k F ) 2 , p d = y (M d ) 2 + (k F ) 2 ; k F) k d F are the Fermi momentums of 
u, d quark, respectively. 

Consequently we can obtain the u, d quark numbers 

rk F pMREk 2 dk 



N u = Vn u = 2N c V 



pMREk^dk 

Jo 2tt 2 ' 1 ' 



N d = Vn d = 2N C V ( F BM^k , (15) 



where V = 4nR 3 /3 is the volume of the droplet, iif = — is the quark number 
density. For given N u and N d , the Fermi momentum of the u, d quark can be determined, 
respectively. 

If the system (droplet) is stable, the thermodynamic potential must satisfy the sta- 
tionary conditions 

on on 

0, 7^ = 0- (16) 



dM u ' dM c 
After some derivation we obtain the constituent quark mass 

f A k 2 E ku dk n „„ . f k F k 2 dk .„„. 
M u = m + 2GN c d m / p MRE k ' u + 2GN c p u d m / Purest , (17) 
Jkf, 2ir z Jo 2ir z 

r A p TP ,AU r k% k 2 clk 

M d = m + 2GN c d m j kd p MRE +2GN c p d d m J o Pmre^ , (18) 

Then we can fix the chemical potentials 

Pu = ^{M u ? + {kp) 2 , p d = ^J(M d ) 2 + (k F ) 2 , (19) 
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and the isospin chemical potential 

A*/ = A* u - A*d • (20) 

By varying the radius i? of the droplet, we can investigate the finite size effect on the 
isospin symmetry breaking and chiral symmetry restoration. 

As a numerical example we take a charge neutral droplet with total quark number 
9000, and the ratio of u, d quark number 1:2. For the current quark mass, we take m = 
6 MeV. For the cutoff and the coupling constant, we take A = 590 MeV and GA 2 = 4.7, 
with which the pion decay constant and pion mass (f n = 93 MeV, = 140 MeV) are 
reproduced well [11J. 
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Figure 1: The dynamical quark masses of u, d quarks as a function of the droplet radius 
R. 



Since the dynamical quark mass can usually be taken to identify the chiral symmetry 
breaking and isospin asymmetry, we investigate the variation of the masses of u and d 
quarks at first. By solving Eqs. ()14|) - (jl8)) . we obtain the dynamical quark mass Mf (f = 
u, d) as a function of radius R of the droplet. The result is illustrated in Fig. 1. From 
Fig. 1 we can notice that the constituent quark mass M u is larger than Md at the same 
radius if the radius of the droplet is neither very small nor extremely large. It means that 
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the isospin symmetry is broken. Whereas with the increase of R, the difference between 
M u and M d becomes small. It indicates that, when the radius R of the droplet is large, the 
quark number density decreases rapidly, and hence the isospin asymmetry becomes small. 
Meanwhile, as the radius decreases (i.e., the quark number density increases since the 
quark number is fixed), the dynamical mass of the quarks can descend to current quark 
mass suddenly. It shows that the chiral symmetry is restored except for the explicit 
breaking. The figure also shows that the critical radius for the chiral symmetry of u 
quark to be restored (R — 17 fm) is smaller than that of d quark (R = 20.5 fm). When 
R is smaller than 17 fm, the dynamical masses of both u quark and d quark descend 
to current quark masses. It manifests that the isospin symmetry can be restored. In 
order to understand the phase transition in a more clear way, we also study the quark 
condensates as a function of baryon number density. The numerical result is illustrated in 
Fig. 2. Fig. 2 shows evidently that the critical density of the chiral symmetry restoration 
of u quark is larger than that of d quark. More concretely, the critical density of d quark 
is smaller than the normal nuclear matter density, whereas, that of u quark is near the 
normal nuclear matter density. All the results show that the contributions of the surface 
and the curvature play important roles for the droplet with a small radius. 
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Figure 2: Quark condensate as a function of baryon number density. (qq) VAC refers to 
the quark condensate with R — > oo. 



7 



40 



> 

CD 

I, 30- 
o 

'■4— ' 

CD 

■g 20- 

Q_ 

ro 
o 



** 
* ★ 

* 



CD 10- 



Q_ 

o 
en 



***************** 



10 



20 
R[fm] 



30 



Figure 3: The isospin chemical potential as a function of the droplet radius R. 



Then, we discuss the isospin chemical potential /i/ in the droplet. By solving Eqs. (14)- 
(20), we obtain the function /// versus the droplet radius. The result is shown in Fig. 3. 
It clearly shows that when R < 17 fm, the isospin chemical potential fij has the minimal 
value 0, and it increases to the maximal value when R = 20.5 fm, then decreases with the 
increase of R. When the radius R is larger than 500 fm, the difference between \i u and is 
less than 1 percentage of the current quark mass. This means that the isospin symmetry 
can be restored when the quark number density is small. The zero isospin chemical 
potential at very small radius shows that the isospin symmetry can also be restored if 
the quark number density is quite large. Another thing important we recognize is that, 
if the droplet takes a size so that only the chiral symmetry of u quark could be restored, 
the isospin chemical potential takes the maximal value about 29.7 MeV in the present 
example of the droplet. It is obvious that such a value is much smaller than the critical 
isospin chemical potential for pion condensate (/if = m n [HI 1211) t° occur. It indicates 
that it may be difficult for pion condensate to happen in the droplet at zero temperature. 

The quark droplets may possibly be formed in the QCD phase transition in the early 
universe or in astro-objects, for example, the main ingredient in the inner part of pulsars 
is quark matter composed of u, d quarks due to the high pressure [213 • Furthermore, 
some works once proposed that the new phase of quark matter may emerge as droplets 
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in nuclear matter at low density [26], even appears as u-d quark stars (P-stars) ^Tf\ l2~Hj. 
Our present results indicate that, as the density of the u-d quark matter in the droplet 
is not very small, the chiral symmetry can be restored, i.e., the chiral phase transition 
takes place. Since the phase transition influences the equation of state of the matter 
drastically, the pressure, the moment of inertia and other characteristics of the matter 
may change suddenly. As shown in Figs. 1 and 2, when the density of the matter in our 
present example is about 0.85n (the radius of the droplet is about 17 fm), the dynamical 
mass of u quark changes from 127 MeV to 6 MeV. Such a sudden change may induce 
a quake and a glitch for the astro-objects. And the light emission strength may also be 
changed abruptly. When the density n G (0.48n , 0.85n ) (corresponding to the present 
example with radius 17 fm< R < 20.5 fm), only the dynamical mass of u quark changes. 
With a further decrease of the baryon number density, the dynamical mass of d quark 
also changes suddenly. Then another glitch may appear in pulsars. These phenomena 
may be taken to identify the effect of chiral phase transition in the droplet of u and d 
quarks |2*9*1 128j. The brightest giant flare from soft gamma-ray repeaters and anomalous 
X-ray pulsars and other observations may also be taken as the observable evidences to 
signal the chiral phase transition effect. It is certain that further studies are necessary to 
carry out a practical analysis. On the other hand, the Figs. 1 and 2 manifest that, in a 
quite large region of nuclear matter density, the isospin symmetry is broken, and hence 
the dynamical mass of u quark is larger than that of d quark. In the view point of bag 
models and soliton models of a hadron, the effective mass of a proton (composing of 2 u 
quarks and 1 d quark) in nuclear matter m* is larger than that of a neutron (consisting 
of 1 u quark and 2 d quarks) m*. Together with the result given in QCD sum rules j50], 
we can recognize that the result m* > m* given in relativistic approaches [HH 1221 EH] has 
solid microscopic foundation. Our present result of the isospin asymmetry provides then 
a clue to solve the controversial problem of the nucleon effective mass splitting in nuclear 
matter [3TJ |321 E3 OH E3 EH] on the theoretical side. It is now rather important since 
experiments have not yet given any conclusion [36]. 

In summary, by taking the two-flavor NJL model and the multiple reflection expansion 
approximation, we have studied the phase transition of a finite size droplet with u and 
d quarks. We find that, the isospin symmetry is preserved and the chiral symmetry is 
broken if the radius of the droplet is extremely large (or the matter density is very small) . 
With the radius changes from infinite to about several hundred fms, both the chiral 
symmetry and the isospin symmetry are broken. In such a case, the dynamical quark 
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mass M u is larger than Md at the same radius R. It provides a clue that the effective 
mass of proton in nuclear matter may be larger than that of neutron. If the radius 
of the droplet is small enough (less than 20.5 fm for N = 9000), the chiral symmetry 
can be restored. And the critical radius for rt-quark is smaller than that for (i-quark. 
Meanwhile, possible observations to identify the chiral phase transition are proposed in 
some astronomical phenomena. In addition, the maximal isospin chemical potential in 
the case of zero temperature is not large enough so as to be comparable with the critical 
isospin chemical potential for pion condensate to emerge. However in our present study, 
we have not taken the flavor mixing J2| into account, where the dynamical quark mass 
Mf (/ = u,d) depends not only on the condensate <pf but also on <pfi. Meanwhile the 
finite temperature effect has not yet been included either. Moreover the experimental or 
astronomical observables identifying the phase transition need also detailed investigations. 
The related studies are under progress. 
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